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Vortex dynamics in molybdenum-germanium superconducting films were found to well approxi- 
mate the unpinned free limit even at low driving forces. This provided an opportunity to empirically 
establish the intrinsic character of free flux flow and to test in detail the validity of theories for this 
regime beyond the Bardeen- Stephen approximation. Our observations are in good agreement with 
the mean-field result of time dependent Ginzburg-Landau theory. 
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I. INTRODUCTION AND BACKGROUND 

The motion of magnetic flux vortices in the mixed state 
of type II superconductors is one of the most studied as- 
pects of superconductors. The most fundamental trans- 
port regime of the mixed state is the state of free flux 
flow (FFF) where the "Lorentz" force on the vortices 
is balanced only by the intrinsic viscous drag, without 
any additional interactions arising from pinning, elastic 
strains, thermal gradients, etc. 

One of the rudimentary characteristics expected of the 
FFF regime (as long as j is not large enough to alter the 
superconducting state) is that the transport response be 
Ohmic, so that the flux-flow chordal resistivity pf = E/ j, 
not just the differential resistivity pd — dE/dj, is con- 
stant; here E is the electric field and j is the current 
density. At sufficiently low values of B, where each vor- 
tex behaves independently, pf is simply proportional to 
the number density of vortices and hence to the mag- 
netic field B. In this limit, elementary theories of flux 
flowi"— have shown that the proportionality factor is of 
the order of Pn/Hc2 resulting in the Bardcen-Stephcn 
(BS) relation pf ^ PnB/Hc2, where p„ is the normal- 
state resistivity and Hc2 is the upper critical field (we 
set ^0 = 1 and use units of tesla for both the magnetic 
field B as well as the magnetizing field H). As B is in- 
creased beyond this low-field limit, additional effects set 
in — such as the suppression of the order parameter with 
increasing B, changes in the circulating current patterns 
around the vortices, et cetera — causing the response to 
deviate from pf (x B. More advanced theories such as 
the work by Larkin and Ovchinnikov^ using a microscopic 
Green-function approach, and various works^^— based on 
time dependent Ginzburg-Landau theory, calculate the 
detailed pf{B,T) behavior beyond the Bardeen-Stephen 
approximation. These theories have not been sufficiently 
verified before. 

While there is an impressive body of experimental work 
that has discovered or confirmed many interesting and 



exotic regimes of vortex dynamics (such as vortex glass, 
melting, flux creep, vortex instability, etc.) the funda- 
mental FFF regime is the least completely studied of 
these. Even the expected Ohmic behavior is usually not 
observed with high accuracy; many previous reportsii"— 
either plot the differential (dE/dj) instead of the chordal 
(E/j) resistivity, and/or use logarithmic scales so that 
the deviation from Ohmic behavior is less conspicuous 
than it can be on a linear-linear plot of p versus j. In 
some earlier workii"— , high current densities were able 
to largely overcome pinning so that p became at least 
roughly constant. These previous works were able to es- 
tablish, at least on a coarse scale, the adherence to the 
basic BS behavior. There are a few partial reportsi^"— 
that compare an observed pf{B,T) curve with the LO 
theory, but the FFF regime has not been systematically 
investigated so as to make a detailed evaluation of the 
different theories that apply to this regime. 

In the present work we were able to cleanly observe 
FFF behavior in low-pinning molybdenum-germanium 
(MoGe) superconducting films, over the temperature 
range 0.6Tc ^ T < Tc, with exact Ohmic behavior visible 
on an uncompressed linear-linear scale even for low cur- 
rent densities {j/jd ^ 1 x 10""' and j/j* r^ Ix 10~^ where 
id and j* are the depairing and flux- flow- instability val- 
ues respectively) that avoid non-linear alterations that 
can arise from high j. Our data allow us to empirically 
elucidate the character of the FFF regime and to test 
theories that make detailed predictions beyond Bardeen- 
Stephen. 

All theories for free flux flow are restricted in one way 
or another (e.g., T close to Tc, B close to Hc2, gap- 
less case, etc.). A microscopic treatment of this phe- 
nomenon, which is unrestricted in its B/Hc2 range (al- 
though restricted to T close to Tc) is that due to Larkin 
and Ovchinnikov- (LO) who flnd (their Eq. 22): 
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where non-linearities that oceur at higher values of B are 
contained in the f{B/Hc2) function that they describe in 
their paper (their Table 1 and related text). 

Alternatively, various authors^"— have calculated the 
flux-flow conductivity using a TDGL (time dependent 
Ginzburg-Landau) approach. The applicable mean-field 
result can be written as: 
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where i^ ~ 0.4 is a roughly T, B, and material inde- 
pendent parameter—. Despite TDGL's strict validity for 
only gapless superconductivity, we find that Eq.[2]fits our 
observations quite well, taking v as an adjustable param- 
eter that we allow to be determined by the experiment. 



II. EXPERIMENTAL TECHNIQUES 

Films of thickness f = 50 nm were sputtered onto 
silicon substrates with 200 nm thick oxide layers from 
a Moo.79Geo.21 alloy target. The deposition system 
had a base pressure of 2 x 10~^ Torr and the argon 
gas working pressure was maintained at 3 mTorr dur- 
ing the sputtering. The growth rate was 0.15 nm/s. 
The samples were patterned into bridges of length I = 
102 /im and width w ~ 6 fim using photolithography 
and argon ion milling. Measurements were conducted 
on four samples A-D. Their transition temperatures Tc, 
normal-state resistances i?„, upper-critical- field slopes 
H'^2 = d,Hc2/dT\T^, corresponding i?c2(0) values (ob- 
tained using the WHH [Werthamcr, Helfand, and Ho- 
henberg] formalismi^) , and Ginzburg-Landau (GL) pa- 
rameter K = 3.54 X 10^-\/— /9„/ioiJ^2 (from Kcs and 
Tsuei^S) are as follows. Sample A: Tc=5.56 K, i?„=555 
17, iy^2=-3-13 T/K, iJc2(0)==12.0 T, and k=78. Sample 
B: Tc=5.41 K, i?„=555 Q, H'^2^-3.13 T/K, Fc2(0)=11.7 
T, and K==78. Sample C: Tc=5.01 K, i?„=630 ft, H'^^^- 
3.0 T/K, i/c2(0)=10.3 T, and k=77. Sample D: Tc=5.00 
K, i?„=540 n, i/^2=-2-63 T/K, i/c2(0)=9.1 T, and k=67. 

The cryostat was a Gryomech PT405 pulsed-tube 
closed-cycle refrigerator that went down to about 3.2 
K. It was fitted inside a 1.3 tesla GMW 3475-50 water- 
cooled copper electromagnet. Calibrated cernox and hall 
sensors monitored T and B respectively. The triggers 
of all measuring instruments were synchronized with the 
pulsed-tube compressor cycle and instrument measure- 
ment windows were set to 1 power line cycle (16.7 ms 
'^3% of the compressor cycle) which ensured a tempera- 
ture consistency of ~10 mK. The main electrical resistiv- 
ity measurements were made using a standard dc four- 
probe method with an in-housc built dc current source 
and voltages measured with a Keithley model 2000 mul- 
timeter. Some extended IV curves were measured using 
0.005 % duty cycle 20 /iS duration pulses (with in-house 
built pulsed-current source and differential preamplifier, 
and a LeGroy model 9314A digital storage oscilloscope) 
to obtain a broader view of the behavior that includes 



higher currents up to and beyond the vortex instability; 
however, these high-j measurements are not essential for 
the analysis and conclusions of the present work, which 
pertains to the low-driving force regime. All data pre- 
sented here were found to be completely reversible and 
showed no hysteresis with respect to cycling of /, B, and 
T . Further details of the measurement techniques have 
been published in previous review papers2ii2^. 



III. RESULTS AND ANALYSIS 

Fig. [1] shows some examples of mixed-state transport 
responses. Panel (a) shows R vs I curves on a linear- 
linear graph to best illustrate the constancy of R. Panel 
(b) shows V vs I curves: the absence of an intercept 
again emphasizes that the response is Ohmic (homoge- 
neously linear) and not just linear with an offset. Both 
panels (a) and (b) were measured with continuous dc 
currents of values that are low compared with jd (by 
10^"*) and j* (by 10~^). Panel (c) shows an example 
of a global view of the transport response that includes 
the high-j behavior measured using pulsed currents. As 
can be seen, the response becomes progressively more 
non-linear with rising j, eventually leading to an insta- 
bility where the curves bend around and enter a nega- 
tive sloped region: There is no reentrant Ohmic behav- 
ior (e.g., please see Fig. 2 of Ref. i21j or Fig. 5 of Ref. HI) 
that one often sees when the low-j Ohmic region is due to 
thermally assisted flux flow^i (TAFF) rather than FFF. 
TAFF, well known from high-temperature superconduc- 
tors, corresponds to the case of a weakly pinned flux 
liquid in the presence of large thermal activation (when 
the pinning energy Uq ^ ksT). TAFF is characterized 
by an exponentially reduced resistivity (as in the case 
of Ref. Il6| where they see an approximately Ohmic re- 
sponse in the pinned regime at feeble current densities 
[10 A/cm^], albeit with a value that is 8% oi Pf) whereas 
here we find p/ of the same order as the nominal Bardeen- 
Stephen value (indeed we find rather detailed agreement 
with Eq. [5]) . When a resistive transition is characterized 
by TAFF, it is expected to show the Arrhenius depen- 
dence logR ~ —1/T] as we show later our resistive tran- 
sitions do not exhibit this behavior but instead follow 
the prediction of Eq. [51 This absence of TAFF behavior 
is consistent with MoGe's small Ginzburg number— i^ 
of Gi = 10-i3K4r2/if^2(0) - 5 X 10-*^. (In Ref. [H in 
the comparable NbGe system, they observe quasi-Ohmic 
behavior presumably corresponding to TAFF; but their 
current densities are 10^ times smaller than ours.) All 
of the Pf data that are subsequently plotted and ana- 
lyzed, were checked for Ohmic behavior and constancy 
of R as in Fig. [Ua), to ensure that the data lie in the 
FFF regime. 

Fig. [2] shows, for samples A (top row) and B (bottom 
row), the free-flux- flow resistances as measured above 
plotted against the magnetic field. The vertical scale cor- 
responds to a normalized resistance range of < R/Rn ^ 
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FIG. 1: (a) R vs I curves for sample A. The horizontal lines 
are least-squares fits to the data, (b) V vs I curves for sample 
B. The lines are least-squares fits to the data extrapolating to 
the origin, (c) Extended V vs I curves (measured with pulsed 
signals) for sample D that include the vortex instability region. 



0.5. (In a later figure we show resistive transitions, which 
include the R -^ i?„ region.) The different curves corre- 
spond to different temperatures as indicated. At low B 
the response is linear and at higher B it curves upward 
(the simple BS formula corresponds to a linear response 
for all B). The left column (panels (a) and (c)) shows fits 
to the LO theory (Eq.[T]) and the right column (panels (b) 
and (d)) shows'fits to the TDGL theory (Eq.E]). While 
both theories are in principle restricted to work only close 
to Tc, the TDGL result provides an excellent description 
of the observed behavior over the entire range, whereas 
the LO fits have shapes that do not conform to the data 
even for T/Tc > 0.9 (top 3 curves) where they ought to 
work (it appears that the / function in Eq. [T] has ex- 
cessive curvature so that the LO fits can't be much im- 
proved even with a scaling constant). Both theories show 
a greater departure when R -^ R„, a region that is bet- 
ter revealed in plots of the resistive transition shown in 
a later figure. For the TDGL fits (panels (b) and (d) of 
Fig- 12): i^ was taken as a fitting parameter for each curve 
so that the dependence i'{T) could be determined. This 
procedure was repeated for the other two samples. The 
resulting values of I'iT) for all four samples are shown 
in Fig. Elja). The observed h'{T) is relatively fiat with 
respect to temperature with values in the 0.2-0.4 range, 
which are consistent with the estimate of i^ ^ 0.4 from 
theoryi^. 

In the literature, we found one comparable R ver- 
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FIG. 2: (a) and (b) Resistance versus magnetic field data 
(symbols) for sample A at temperatures (bottom to top): 
T=3.25, 3.74, 4-24, 4.56, 4.79, 4.90, 5.00, 5.21, 5.31 K. (c) 
and (d) Resistance versus magnetic field data (symbols) for 
sample B at temperatures (bottom to top): T — 3.17, 3.48, 
3.71, 4.01, 4.21, 4.53, 4.77, 4.88, 4.98, 5.19, 5.30 K. (a) 
and (c) show fits to LO theory (Eq. Qp. (b) and (d) show fits 
to the TDGL theory (Eq. \^ using u as a fitting parameter 
for each curve (i.e., for each T). The resulting values off 
are shown in Fig. [3 All fits use the same Hc2 (T) given by 
Hc2{T) = H'^2lT - Tco] and the measured value H'^2=-3.125 
T/K. 



sus B curve for another low pinning system (amorphous 
NbsGe) by Berghuis et alJ^ (even though they plot dif- 
ferential rather than the chordal resistivity, their data 
appear to lie approximately in the FEE regime). Their 
data are shown in Eig.[3l^b). As can be seen, the solid red 
line corresponding to the TDGL function (with v — 0.27) 
fits the entire range of their data well (this v is shown 
in panel (a) as an asterisk). In their paper, the authors 
fit the data to the LO prediction for the B k, Hc2 region 
(Eq. 30of LO): 



O-f = O-n + 0-„a[l - B/Hc2] 



(3) 



with a = 2.44. That fit is shown as a dashed blue line 
in Eig. ^h) (in their paper the authors show only the 
first-order linear expansion Pf/pn ~ 1 — ^[l — B / Hc2\ of 
Eq.[3]). On the same plot we have also shown a curve cor- 
responding to Eq. [T] (the LO prediction for T fa T^ since 



these data correspond to T/Tc = 0.85) as the black dot- 
ted line. As can be seen, both LO curves show marginal 
agreement with the data and only over a very narrow 
range. 
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FIG. 3: (a) Experimentally deduced parameter v (of Eq. \^ 
and its variation with temperature for MoGe samples A-D. 
The asterisk shows a v value for the Nbs Ge data of refer- 
ence\l^ plotted in panel (b). (h) Normalized resistance versus 
normalized field for amorphous Nbz Ge films from Berghuis et 
alJ^. t = T/Tc = 0.85 with Tc^2.93 K. Solid red line repre- 
sents a TDGL curve (Eq. \^ with v=0.27. Blue dashed and 
black dotted lines correspond to the LO theory for the condi- 
tion "close to Hc2 " (Eq. 0) with a=2.44 cii^d for the condition 
"close to Tc" (Eq.\Vj) respectively. 



To explore the region close to T k, Tc and Rf ~ Rn, we 
show in Fig.|4]the resistive transitions in various magnetic 
fields. Panel (a) shows the data for sample A along with 
LO curves (Eq. [1} ; there is some resemblance in trends 
between the data and theory, but for the most part the 
agreement is not close. Panels (b) and (d) show the data 
for samples A and B along with TDGL curves (Eq. [2|) . 
Given our earlier empirical confirmation that v is T in- 
dependent (Fig. [3l[a)), we tried to fit all curves for each 
sample with a single v value. As can be seen, for each of 
panels (b) and (d), the TDGL theory provides excellent 
agreement with the observed behavior over the lower and 
middle portion of the Rf / Rn range with a single parame- 
ter. The obtained values of i' = 0.34, 0.26, 0.23, and 0.31 
(for samples A-D respectively) are quantitatively consis- 
tent with the theoretically estimated v ^ 0.4. The region 
right around the transition does not fit either theory well, 
perhaps because of a cross over to 2 dimensionalityi^ or 
additional effects^^ not included in the present theories, 
as well as because of the finite width of the transition. 



Panel (c) plots the data for sample A as log R versus 1/T. 
As can be seen it does not exhibit the Arrhenius behavior 
{R ex e^^/^) that characterizes the TAFF regime. 
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FIG. 4: Resistive transitions in magnetic fields. Symbols show 
experimental data. Solid lines are theoretical curves, (a) Sam- 
ple A with LO theory curves (Eq. Qp. (b) Sample A with 
TDGL theory curves (Eq. [1 with u=0.34 for all B and T). 
(c) Arrhenius plots of the same resistive-transition data for 
Sample A. (d) Sample B with TDGL theory curves (Eq. [1| 
with v=0.26 for all B and T). 



IV. SUMMARY AND CONCLUSIONS 

Free fiux fiow is the most primitive regime of transport 
in the mixed state; however, its experimental investiga- 
tion is normally hindered by gross alterations of the dy- 
namics due to pinning and other complications (dissocia- 
tion into pancakes due to a layered structure, exotic gap 
symmetries, melting/entanglement due to high tempera- 
tures, etc.). The present experiment has observed clean 
FEE behavior in one of the simplest and nearly model 
superconductors (unpinned, isotropic, low-temperature, 
weak-coupling BCS, etc.). This allowed us to look closely 
at the transport characteristics of this regime beyond the 
simple Bardecn-Stephen formula pf ^ PnB / Hc2 and to 
assess the applicability and accuracy of the microscopic 
and TDGL theories. We find that the mean-field result 



arising out of TDGL provides a much better agreement 
with the observed pf{B,T) function — and is apphcable 
over a wider range of conditions — than the theory of 
Larkin and Ovchinnikov. 

Besides the one experimental curve from Ref. llSl 
(shown in our Fig. Eljb)) where the authors apply the 
near-_ffc2 LO Pf{B, T) function (Eq.[31 shown as the blue 
dashed line in Fig. [3l^b)), in the works by Babic et al. on 
vortex noisei^ and vortex instabilitj*ii, the authors show 
(incidentally to the main investigations of those papers) 
partial agreement with the near-T,, LO pf{B,T) func- 
tion (Eq.^; the insets of their Figs. 1 of both Ref. [la 
and Ref. |l7| each show one experimental curve for one 
temperature for one sample. 

The present work only evaluates the accuracy of the 
Pf{B, T) function in the FFF regime and makes no com- 
ments on the accuracy and applicability of the LO theory 
for other regimes such as the vortex instability. An im- 
portant value of a microscopic theory such as LO lies in 
its ability to go beyond phenomenology and make a con- 
nection between measurements and microscopic param- 
eters (such as the extraction of the inelastic scattering 
time as carried out in Ref. KLSl ) . 

On the other hand, the mean-field result of the TDGL 



theory provides a rather excellent account of the observed 
Pf{B,T) function as seen from Figs. [2] and 2] along with 
quantitative agreement with the predicted parameter v. 
This is a useful observation, because in some calcula- 
tions of exotic effects — such as vortex instabilities — the 
BS formula is often used as a starting point. For such ap- 
plications, the TDGL formula (Eq. ^ with i/ - 0.2-0.4 
may provide a more accurate alternative while retaining 
most of the simplicity of the BS result. 
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